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MAGIC J0616+225 as delayed TeV emission of cosmic-rays 
diffusing from SNR IC 443 
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ABSTRACT 

We present a theoretical model that explains the high energy phenomenology of the neighbor- 
hood of SNR IC 443, as observed with the Major Atmospheric Gamma Imaging Cherenkov 
(MAGIC) telescope and the Energetic Gamma-Ray Experiment Telescope (EGRET). We in- 
terpret MAGIC 10616-1-225 as delayed TeV emission of cosmic-rays diffusing from IC 443 
and interacting with a known cloud located at a distance of about 20 pc in the foreground of 
the remnant. This scenario naturally explains the displacement between EGRET and MAGIC 
sources, their fluxes, and their spectra. We compare this model with others recently presented, 
and discuss how it can be tested with observations by the Gamma-ray Large Area Telescope 
(GLAST). 
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1 INTRODUCTION 

Recently, MAGIC presented the results of observations towards 
SNR IC 443, yielding to the detection of a new source of y- 
rays, J0616+225 (Albert et al. 2007). This source is located at 
(RA,DEC)=(0646'"43% +22°31'48"), with a statistical positional 
error of 1.5', and a systematic error of 1'. A simple power law 
was fitted to the measured spectral points: dA'y/(dAdfd£') = (1.0 ± 
0.2) X 10-" (£'/0.4TeV)"^ '*''^cm-2s-'TeV-', with quoted errors 
being statistical. The systematic error was estimated to be 35% in 
flux and 0.2 in spectral index. No variability was found along the 
observation time (over one year). No significant tails nor extended 
structure was found at the MAGIC angular resolution. These results 
were confirmed by observations with the VERITAS array (Humen- 
sky et al. 2008). In addition, consistent upper limits were reported 
by Whipple (Holder et al. 2005) and CAT (Khelifi 2003). 

MAGIC J0616+225 is displaced with respect to the posi- 
tion of the non- variable (Torres et al. 2001) EGRET source 3EG 
I0617-I-2238 (Hartman et al. 1999). Indeed, the EGRET central po- 
sition is located directly towards the SNR, whereas the MAGIC 
source is south of it, close to the 95% CL contour of the EGRET 
detection. As Albert et al. (2007) showed, the MAGIC source 
is located at the position of a giant cloud in front of the SNR, 
it would not be surprising if they are related, which we explore 



here. The EGRET flux is (51.4±3.5) xlO » ph cm" 



with 



a photon spectral index of 2.01±0.06. An independent analy- 
sis of GeV photons measured by EGRET resulted in the source 
GeV J0617+2237 (Lamb & Macomb 1997), also at the same lo- 
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cation of 3EG I0617-I-2238 at the center of the SNR. Extrapolating 
the spectrum of the EGRET source into the VHE regime, we would 
obtain a higher flux and a harder spectrum than which was observed 
for MAGIC 10616+225, supporting the view that a direct extrapo- 
lation of this and other EGRET measurements into the VHE range 
is not valid (Funk et al. 2008). 

Here we present a theoretical model (see Aharonian & Atoyan 
1996, also Gabici & Aharonian 2007) explaining the high en- 
ergy phenomenology of IC 443, focussing in the displacement be- 
tween EGRET and MAGIC sources. Our interpretation of MAGIC 
10616+225 is that it is delayed TeV emission of cosmic-rays (CRs) 
diffusing from the SNR. We compare this model with others re- 
cently presented, and discuss how it can be tested with observations 
with the Gamma-ray Large Area Telescope (GLAST). 



2 THE SNR IC 443 IN CONTEXT OF THE MODEL 

IC 443 is an asymmetric shell-type SNR with a diameter of ~45 
arc minutes (e.g., Fesen & Kirshner 1980). Two half shells appear 
in optical and radio images (e.g, Braun & Strom 1986, Leahy 2004, 
Lasker et al. 1990). The interaction region, with evidence for mul- 
tiple dense clumps, is also seen in 2MASS images (e.g. Rho et al. 
2001). In radio, IC 443 has a spectral index of 0.36, and a flux 
density of 160 Jy at 1 GHz (Green 2004). Claussen et al. (1997) 
reported the presence of maser emission at 1720 MHz at (/, b) ~ 
(-171.0, 2.9). Recently, Hewitt et al. (2006) confirmed Claussen's 
et al. measurements and discovered weaker maser sources in the 
region of interaction. IC 443 is a prominent X-ray source, observed 
with Rosat (Asaoka & Aschenbach 1994), ASCA (Keohane 1997), 
XMM (Bocchino & Bykov 2000, 2001, 2003, Bykov et al. 2005, 
Troja 2006), and Chandra (Olbert et al. 2001, Gaensler et al. 2006). 
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The works by Troja et al (2006) and Bykov et al. (2008) summa- 
rize these observations. In what follows we present some additional 
features of IC 443, relevant for our model. 

Age: A small (~ 1000 yr) age was put forward by Wang et al. 
(1992) and seconded by Asaoka & Aschenbach (1994) and Keo- 
hane (1997), although IC 443 is now agreed to have a middle- 
age of about 3 X 10'' yrs. This age has been initially advocated by 
Lozinkaya (1981) and was later consistently obtained as a result 
of the SNR evolution model (Chevalier 1999). Recent observations 
by Bykov et al. (2008) confirm that there are a few X-ray-emitting 
ejecta fragments, a number much smaller than which should be the 
case in a younger SNR. 

Distance: Kinematical distances from optical systemic veloc- 
ities span from 0.7 to 1.5 kpc (e.g., Lozinskaya 1981). The assump- 
tion that the SNR is associated with a nearby HII region, S249, im- 
plies a distance of ~ 1.5 - 2.0 kpc. Several authors claimed that the 
photometric distance is more reliable (e.g, Rosado et al. 2007), and 
concurrently with all other works on IC 443, we adopt its distance 
as 1.5 kpc (thus, 1 arcmin corresponds to 0.44 pc). 

Energy of the explosion: There is no clear indicator for E^i, 
the energy of the explosion in units of 10^' erg. Chevalier (1999) 
obtains a lower limit of 4 x 10'" erg, whereas lower estimations are 
provided by Dickman (1992), based on Mufson (1986), -albeit the 
latter assumed an age of -5000 yr. Lacking a strong reason for other 
numerical assumption, we will assume that £51 = 1, although to 
be conservative, we will subsequently assume that only 5% of this 
energy is converted into relativistic CRs. Reasonable differences in 
our assumed value of £51 are not expected to have any impact on 
this model. 

The molecular environment: Comett et al. (1977) and De 
Noyer (1981) were among the first to present detailed observa- 
tions of molecular lines towards IC 443. Subsequently, Dickman 
et al. (1992), Seta et al. (1998), Butt et al. (2003) and Torres et al. 
(2003) among others, presented further analysis. These works con- 
form the current picture for the environment of IC 443: a total mass 
of ~ 1.1 X 10"* Mq mainly located in a quiescent cloud in front of 
the remnant (with linear scales of a few parsecs and densities of a 
few hundred particles cm"') that is absorbing optical and X-ray ra- 
diation (e.g., Lasker 1990, Troja et al. 2006), a scenario aheady put 
forward by Comett et al. (1977). Dickman et al. (1992) estimated 
that 500-2000 Mq are directly perturbed by the shock in the north- 
ern region of interaction, near the SNR itself. Huang et al. (1986) 
found several clumps of molecular material along this interacting 
shell, with subparsec linear scales. Rosado et al. (2007) found in- 
homogeneities down to 0.007 pc. As it is usual, we will neglect 
these latter inhomogeneities when considering the propagation of 
CRs in the ISM, i.e. we thus assume an homogeneous medium of 
typical ISM density where CRs ditfuse. Then, the molecular mass 
scenario is a main giant cloud in front of the SNR containing most 
of the quiescent molecular material found in the region, and smaller 
cloud(s) totalizing the remaining mass located closer to the SNR. 



3 DIFFUSION OF CRS FROM IC 443 

The spectrum of y-rays generated through ;;r''-decay at a source 
of proton density is Fy{Ey) = 2 f^^^JF^iE;,)/ ^jEl- ml) dE„, 

where £J"'"(£'^) = Ey + mllAEy, and F„{E„) 

^^^P CL Jp{E){dcr„{E„, Ep)/dE,} dE^. Here, da^E^, E,)ldE, 

p 

is the differential cross-section for the production of ^i^-mesons of 
energy E^ by a proton of energy Ep in a pp collision. For an study 



of different parameterizations of this cross section see Domingo- 
Santamaria & Torres (2005) and Kelner et al. (2006). The limits 
of integration in the last expression are obtained by kinematic 
considerations (see e.g., Torres 2004). In these expressions we have 
implicitly assumed a uniform cosmic-ray density in the cloud as 
well as in the cloud's gas number density (the size of the molecular 
cloud is smaller than the distance to IC 433, we therefore neglect 
the temporal, spatials effects within the molecular cloud itself; the 
whole molecular clouds becomes instantly a cosmic ray target). 

The CR spectrum is given by Jp(E, r, t) = [cfi/An]/, where 
f{E, r, t) is the distribution function of protons at an instant t 
and distance r from the source. The distribution function satisfies 
the radial-temporal-energy dependent diffusion equation (Ginzburg 
& Syrovatskii 1964): (df/dt) = {D(E)/r^)(d/dry(df/dr) + 
(d/dE) (Pf) + Q, where P = -dE/dt is the energy loss rate of the 
particles, Q = Q(E, r, t) is the source function, and D{E) is the dif- 
fusion coefficient, for which we assume here that it depends only on 
the particle's energy. The energy loss rate are due to ionization and 
nuclear interactions, with the latter dominating over the former for 
energies larger than 1 GeV. The nuclear loss rate is Prmc = E/Tpp, 
with Tpp = (rip CKCTppY^ being the timescale for the corresponding 
nuclear loss, k ~ 0.45 being the inelasticity of the interaction, and 
CTpp being the cross section (Gaisser 1990). Aharonian & Atoyan 
(1996) presented a solution for the diffusion equation for an arbi- 
trary energy loss term, diffusion coefficient, and impulsive injection 
spectrum /inj(-E), such that Q(E, r, t) = Nofmj(E)Sr5(t). For the par- 
ticular case in which D(E) oc E^ and /i„j oc E'", the general solution 
is /(£, r, t) ~ (NoE-VTt^'^Rl^) exp [-{a - l)t/Tpp - (R/R^f] , 
where R^.f = 2(D(E)t[exp(tS/Tpp) - l]/[tS/Tpp]y'^ stands for the 
radius of the sphere up to which the particles of energy E have 
time to propagate after their injection. In case of continuous in- 
jection of accelerated particles, given by Q(E, t) = QoE'''T{t), 
the previous solution needs to be convolved with the function 
T(t - t') in the time interval ^ < f. If the source is de- 
scribed by a Heavside function, T(t) = 0(0 Atoyan et al. (1995) 
have found a general solution for the diffusion equation with ar- 
bitrary injection spectrum, which with the Usted assumptions and 
for times t less than the energy loss time, leads to: f{E, r, t) = 
(QoE-" /47TD(E)r)(2/ y/^) H e-'^dx. We will assume that a = 
2.2 and make use of these solutions in what follows. 

Fig. 1 shows the current CR spectrum generated by IC 443 
at two dilferent distances from the accelerator, 10 (solid) and 30 
(dashed) pc. The SNR is considered both as a continuous acceler- 
ator with a relativistic proton power of Lp = 5 x 10^' erg s~' (the 
proton luminosity is such that the energy injected into relativistic 
CRs through the SNR age is 5 x lO'*' erg), and an impulsive in- 
jector with the same total power (injection of high energy particles 
occur in a much shorter time than the SNR age). The horizontal line 
in Fig. 1 marks the CR spectrum near Earth, so that the excess of 
CRs in the SNR environment can be seen. For this example, the dif- 
fusion coefficient at 10 GeV, Dio, was chosen as 10'^ cm^ s"', with 
S = 0.5. CRs propagate through the ISM, assumed to have a typical 
density. In the scale of Fig. I, curves for riisM = 0.5, 1,5, and 10 
cm"' would be superimposed, so that njsM becomes an irrelevant 
parameter in this range (this stems from the fact that the timescale 
for nuclear loss Tpp obtained with the densities considered for the 
interestellar medium, nisM, is orders of magnitude larger than the 
age of the accelerator). Differences between the different kind of 
accelerators assumed are also minimal for the SNR parameters. 

Fig. 2 shows the result for the y-ray emission coming from 
the cloud located at the position of the MAGIC source, when we as- 
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Figure 1. Cuirent CR spectrum generated by IC 443 at two different dis- 
tances, 10 (solid) and 30 (dashed) pc, at the age of the SNR. Two types of 
accelerator are considered, one providing a continuous injection (black) and 
other providing a more impulsive injection of CRs (red). The horizontal line 
marks the CR spectrum near the Earth. The y-axis units have been chosen 
to emphasize the excess of CRs in the SNR environment. 



sume it lies at different distances in front of IC 443. The giant cloud 
mass is assumed (consistently with observations) as 8000 Mq. The 
accelerator properties and power of IC 443 are as in Fig. 1, in each 
case. Fluxes are given for an ISM propagation in a medium of ;i = 1 
cm"', although again we have checked this is not a relevant param- 
eter as discussed above. We find that clouds located from ~20 to 
~30 pc produce an acceptable match to MAGIC data. In the case of 
a more impulsive accelerator, the VHE predicted spectra is slightly 
steeper than that produced in the continuous case at the same dis- 
tance, so that it provides a correspondingly better fit to the MAGIC 
spectrum. Fig. 2 also shows, apart from MAGIC data, EGRET mea- 
surements of the neighborhood of IC 443. We recall that these two 
sources are not located at the same place, what we emphasize us- 
ing different symbols. Fig. 2 shows that there is plenty of room 
for a cloud the size of that detected in front of IC 443 to gener- 
ate the MAGIC source and not a co-spatial EGRET detection. In 
the case of GLAST, measurement of this region will allow us to 
constrain the separation between the SNR and the cloud, since for 
some distances a GLAST detection is also predicted. The existence 
of a VHE source without counterpart at lower energies is the result 
of diffusion of the high-energy CRs from the SNR shock, which is 
an energy dependent process leading to an increasing deficit of low 
energy protons the farther is the distance from the accelerator. 

To clarify our previous assertion, and since our solution to the 
diffusion-loss equation is a function of time, we show the evolution 
of the flux along the age of the SNR. In Fig. 3 we show the inte- 
grated photon flux coming from the position of the giant cloud as 
a function of time above 100 MeV and 100 GeV in the impulsive 
case. Different qualities of the accelerator (impulsive or continu- 
ous) produce a rather comparable picture. At the age of the SNR 
(the time at which we observe) GLAST should see a source only 
for the closest separations. On the contrary, the integrated photon 
fluxes above 100 GeV present minimal deviations, and a MAGIC 
source is always expected. 

Fig. 2 also presents the results of our theoretical model focus- 
ing in the energy range of EGRET. There, the CR spectrum interact- 
ing with a local-to-the-SNR cloud is obtained assuming an average 
distance of interaction of 3^ pc. A few hundred Mq located at this 
distance (~700 Mg for the case of an impulsive, and ~300 Mq for 
a continuous case) produce an excellent match to the EGRET data. 
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Figure 2. MAGIC and EGRET measurement of the neighborhood of 
IC 443 (stars and squares, respectively) as compared with model pre- 
dictions. The top (bottom) panel shows the results for an impulsive 
(continuous) case. At the MAGIC energy range, the top panel curves 
show the predictions for a cloud of 8000 Mq located at 20 (1), 25 (2), 
and 30 (3) pc, whereas they correspond to 15 (1), 20 (2), 25 (3), and 
30 (4) pc in the bottom panel. At the EGRET energy range, the curve 
shows the prediction for a few hundred M© located at 3^ pc. The 
EGRET sensitivity curve (in red) is shown for the whole lifetime of the 
mission for a typical position in the Inner Galaxy, dominated by diffuse 
y-ray background. The GLAST sensitivity curve (in blue) (taken from 
http://www-glast.slac.stanford.edu/software/IS/glast_latperformance.html I 
show the 1-year sky-survey sensitivity for the same position in the Iimer 
Galaxy. 

without generating a co-spatial MAGIC source. Concurrently with, 
e.g., Gaiser et al. (1999), we find that the lowest energy data points 
in the EGRET range are produced by bremsstrahlung of accelerated 
electrons, curves that for simplicity we do not show in Fig. 2. 

As spinoff of the constraints provided by the observed phe- 
nomenology (e.g., the molecular environment and the position of 
the y-ray sources) in the setting of this model, we find that Dio 
should be low, of the order of 10^'' cm^ s"' . By varying the diffusion 
coefficient and studying its influence in our results, we obtain that 
if the separation between the giant cloud and the SNR is >10 pc, 
an slower diffusion would not allow sufficient high energy particles 
to reach the target material; thus, the MAGIC source would not be 
there. On the other hand, if the separation between the main cloud 
and the SNR is <10 pc, we would have detected an EGRET source 
at the position of the cloud, which is not the case. Our model, there- 
fore, allows to put constraints on Ojo at 1.5 kpc from Earth, near 
IC 433, using the fact that the MAGIC and EGRET sources are dis- 
placed from each other. Such values of Dio are expected in dense 
regions of ISM such as the one we study (Ormes et al. 1988, Gabici 
& Aharonian 2007). 
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Figure 3. Integrated photon flux as a function of time above 100 MeV and 
100 GeV, solid (dashed) lines con'espond to the case of the cloud located at 
10 (30) pc from an impulsive accelerator. The horizontal lines represent the 
values of integrated fluxes in the case that the CR spectrum interacting with 
the cloud is the one found near Earth. The vertical line stands for the SNR 
age. EGRET and GLAST integral sensitivity, consistent in value and color 
coding with those in Fig. 2, are shown. 



4 DISCUSSION 

Here we put our model in the context of others. Bykov et al. (2000) 
suggested that the GeV emission seen towards IC 443 is mostly 
due to relativistic bremsstrahlung. As noted by Butt et al. (2003), 
who already favor a hadronic emission of y-rays, the synchrotron 
radio emission seen towards the rim of the SNR and the centrally 
located EGRET source must then be associated. Judging from the 
localization of the multiwavelength emissions, this does not seem 
to be the case. An unavoidable (but subdominant at high energies) 
bremsstrahlung component coming from primary and secondary 
electrons can however play a role at the lowest energies in the 
EGRET range, as already shown by Gaisser et al. (1999) and found 
above. 

Bocchino & Bykov (2001) suggested that the systematic er- 
rors in the EGRET location contours could yield the pulsar CXOU 
J061705. 3+222127 discovered by Olbert et al. (2001) and its neb- 
ula, as the source of the GeV emission. As discussed above, this 
contradicts current observations. 

The same assumption (i.e., that the positions of the measured 
EGRET and GeV sources are wrong by half a degree) was taken 
by Bartko & Bednarek (2008). The pulsar nebula is also displaced 
from the MAGIC detection by 20 arcmin, but these authors sug- 
gested that they may be connected if a pulsar with a velocity of 
250 km s"' moves along the SNR age, i.e., the pulsar CXOU 
J061705. 3-1- 222127 should have been borne at the SNR center and 
travel to its current position while accelerating particles that inter- 
act with the cloud, giving rise to the MAGIC source. We do not 
see that this scenario matches the observed phenomenology: the 
EGRET source should be on top of the current position of the pul- 
sar and not where it actually is, and physically, it should be the 
result of pulsed emission (like in Vela), although pulses were not 
reported. The only argument supporting the latter assumption is 
that the flux and spectrum of 3EG J0617-I-2238 is similar to that 
of PSR 1706-44, observed by EGRET. This would apply to dozens 
of other EGRET sources and can not be sustained as circumstan- 
tial evidence of physical similarity (e.g., Torres et al. 2001b, 2003, 
Reimer 2001, Romero et al. 1999). In addition, the MAGIC source 
is generated by inverse Compton from electrons accelerated at an 
initial phase of the pulsar and traveling towards the cloud. Since 



the difference in target photon fields in the region surrounding the 
cloud should not be significant, and the target field should even be 
larger at the position of the interacting shock in the northeast, the 
localization and size of the MAGIC source is not explained. 

Very recently, Zhang and Fang (2008) presented an alternative 
model for IC 443, in which a fraction of the SNR shell evolves in 
the molecular cloud, and other in the ambient interstellar environ- 
ment, encountering different matter densities. In this model, the y- 
rays observed by EGRET are mainly produced via pp interactions 
with the ambient matter in the clouds, as is the MAGIC source. 
Although this may sound similar to our model, the former is a key 
difference between them: in Zhang and Fang scenario, both EGRET 
and MAGIC sources should be at the same position. This is a con- 
sequence of the fact that in this model, the radial dependence of 
the CR spectrum is not considered. As Gabici & Aharonian (2007) 
noted in general, and as we have seen in our IC 443 results above, 
an old SNR cannot confine multi-TeV particles in their shells. 



5 CONCLUDING REMARKS 

Here we have shown that MAGIC J0616-I-225 is consistent with 
the interpretation of CR interactions with a giant molecular cloud 
lying in front of the remnant, producing no counterpart at lower 
energies. We have also shown that the nearby EGRET source can 
be produced by the same accelerator, and that in this case, a co- 
spatial MAGIC source is not expected. In our model, the displace- 
ment between EGRET and MAGIC sources has a physical origin. 
It is generated by the different properties of the proton spectrum at 
different locations, in turn produced by the diffusion of CRs from 
the accelerator (IC 443) to the target. Specific predictions for fu- 
ture observations can be made as a result of this model. At high 
energies, we should see a morphological and spectral change from 
the position of the cloud (i.e. the center of MAGIC J0616-H225) 
towards the center of IC 443. At a morphological level, the lower 
the energy, the more coincident with the SNR the radiation will be 
detected. At a spectral level: sufficient statistics should show that 
the lower the y-ray energy the harder the spectrum is. Both pre- 
dictions should show in future MAGIC II observations, and as a 
combination of GLAST and MAGIC data. GLAST observations, 
in addition, may be sensitive enough to detect the same cloud that 
shines at higher energy, which ultimately will allow to determine its 
separation from the remnant, if the diffusion coefficient is assumed 
-as we showed-, or viceversa. 
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